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Localization problem of electronic states in a two-dimensional quantum spin Hall system (QSH - a 
symplectic model with a non-trivial topological structure) is studied by the transfer matrix method. 
The phase diagram in the plane of energy and disorder strength is exposed, and demonstrates 
"levitation" and "pair-annihilation" of the domains of extended states analogous to that of the 
integer quantum Hall system. The critical exponent v for the divergence of the localization length is 
estimated as v = 1.6 which is distinct from both exponents pertaining to the conventional symplectic 
and the unitary quantum Hall systems. This strongly suggests a different universality class related 
to the non-trivial topology of the QSH system. 
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PACS numbers: 73.20.Fz, 72.25.-b, 73.43.-f, 03.65.Vf 

The spin Hall effect (SHE) in semiconductors, namely, 
generation of spin current perpendicular to an applied 
electric field, has recently attracted intensive attention 
following its theoretical proposals Recent exper- 

iments observed an accumulated spin polarization near 
the edges of a sample in two-dimensional n-type Q and p- 
type |fjj quantum wells of GaAs under an applied electric 
field. There still remain some controversial issues about 
the origin of these accumulated spins Jo, i.e., whether it 
is due to extrinsic 0, H, or intrinsic [JQ origins. How- 
ever, the fundamental importance of the SHE is the un- 
derlying topological structure of Bloch wavefunctions in 
systems with time-reversal (TR) symmetry. The concept 
of SU(2) non-Abelian Berry phase plays a crucial role in 
the intrinsic SHE [l|, |2j , which is a generalization of the 
anomalous Hall effect in metallic ferromagnets driven by 
the U(l) Berry phase [lfj. These topological aspects led 
to the concept of spin Hall insulator (SHI) which is 
a band insulator showing nonzero spin Hall conductivity. 
In these models edge modes are shown to be gapful 
in generic cases |l2( , in sharp contrast with quantum Hall 
systems (see also ref. where another proposal for the 
quantum spin Hall (QSH) system is discussed). 

Originally, zero gap semiconductors such as HgTe, a- 
Sn and narrow gap semiconductors such as PbTe were 
listed as candidates for SHI 11]. Later on, a model for 
QSH effect in graphene has been proposed [lj,lll|> which 
is shown to be topologically distinct from the previous 
ones since it is characterized by an odd number Z2 of 
pairs of gapless edge modes for a semi-infinite system. In 
that case, Kramers theorem for systems with TR symme- 
try prevents the hybridization of edge modes with oppo- 
site chiralities 0, 0, ^| . The number Z2 so defined is 



claimed to be related to the mod 2 index of the real Dirac 
operator and the latter is then explicitly constructed from 
the Bloch wavefunctions and used to distinguish the QSH 
phase from a simple SHI. Note however that spin is not 
a conserved quantity in the presence of spin-orbit inter- 
action, and one cannot define a topological invariant of 
the first Chern class for the corresponding SU(2) gauge 
field B . (There is one proposal for the conserved spin 
current [18j which applies only for spatially homogeneous 
systems.) Therefore, spin Hall conductance is not related 
to a topological integer such as the U(l) Chern number 
Ctj(i) for the charge quantum Hall effect (QHE) |l9j. In 
this context, a recent work introduces a spin Chern num- 
ber C sc (and a Chern number matrix), and investigates 
its stability toward disorder and its relationship to edge 
modes [13 0, |52]. An important consequence drawn 
therein is that even though there is no conserved spin 
current, one can still define a conserved topological num- 
ber associated with twisted (spin-dependent) boundary 
conditions [2^] . Here the analogy with the QHE is rather 
appealing. 

A question of paramount interest is how the occur- 
rence of these new topological quantities affect the local- 
ization/delocalization properties of wave functions in the 
presence of disorder. Although the standard symplec- 
tic ensemble exhibits a metal-insulator transition in two- 
dimensions, here we encounter a symplectic system with 
a non-trivial topological structure. It is already known 
for the quantum Hall systems, i.e., the unitary case, that 
Cu(i) [Hi associated with the conserved charge current 
protects the isolated extended states at some special ener- 
gies [2^ , although the naive renormalization group (RG) 
study of the unitary universality class dictates the local- 
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ization for any finite disorder strength in two-dimensions 
[25| . Since Cum is defined as a response to twisted 
boundary conditions, existence of extended states is nec- 
essary for having a nonzero Cum ■ This topological inte- 
ger cannot change continuously. Rather, it jumps when 
two extended states with opposite Chern numbers an- 
nihilate each other. The pertinent field theory is the 
nonlinear-er model with a topological term, which results 
in two-parameter scaling RG equations Assuming 
that Z2 classification or C sc is also associated with ex- 
tended states, we expect a similar scenario for the local- 
ization problem in the symplectic model belonging to the 
non-trivial Z2 classification and/or with nonzero C' sc . 

Thus, in this work, we investigate the localization 
problem for a symplectic model with a nontrivial topolog- 
ical structure, i.e., QSH system with odd Z2 and C sc = 2. 
Using the transfer matrix method, the phase diagram 
in the plane of disorder strength and energy is revealed, 
which manifests the features of levitation and pair annihi- 
lation of extended states similar to the unitary QHE case, 
albeit with finite energy width of the extended region. Fi- 
nite size scaling analysis of the localization/delocalization 
transition yields an exponent v = 1.6 for the divergence 
of the localization length, which is distinct from both 
that of the symplectic {y = 2.73) [2(| and that of the 
unitary QHE (y = 2.33) universality classes. This 
strongly suggests a new universality class of the localiza- 
tion/delocalization transition in the symplectic class, due 
to the existence of a non-trivial topological structure. 

The model Hamiltonian we study is a generalization of 
that proposed for graphene by Kane and Mele [Til fl5j | : 



H = 



\t + iX K (a x dij) z 

(H) 

^ c\{t' + i\soVijcr z )cj 
(y)' 

^ cl(X v rii +Wi + h x cr x )c, 



(1) 



Here c"p is the spinor annihilation (creation) operator, <x 
is the set of the Pauli matrices, t is the conventional hop- 
ping energy between nearest neighbor sites (ij) (we use 
t = 1 as the unit of energy), Ar is the Rashba spin-orbit 
interaction strength and is the unit vector connect- 
ing {ij}. The constant t' is the conventional hopping 
energy between second nearest neighbor sites {ij}' ■ Here 
it is introduced just for assuring stability of the numerical 
analysis by the transfer matrix method and is fixed at the 
small value 0.01. Moreover, Aso represents the spin-orbit 
interaction strength with — (2/v3)(di X d-i) = ±1, 
where d\ and d2 are unit vectors along the two bonds 
connecting {ij}' , while X v is the alternation of the site 
energies between the A and B sublattices (r^ = ±1). 
The random potential is uniformly distributed be- 
tween — W/2 and W/2. Finally, the last term represents 



the magnetic field along the x-direction, which breaks TR 
symmetry. The phase diagram of this model without the 
random potential Wi and the magnetic field h x has been 
already displayed in the inset of Fig. 1 in ref. [1J| . In the 
plane Au/Aso-Ar/Aso, there is a finite domain of QSH 
state with Zi = 1, C sc = 2, bounded by a curve where 
the gap closes. Outside this domain, the gap opens up 
again and the system becomes the usual SHI with Z2 = 0, 
C sc = 0. Note that for Ar = 0, the system is decou- 
pled into two independent unitary subsystems H^, |2^| . 
For X v < A£, each unitary model has Cum's of oppo- 
site signs for the valence and conduction bands, and the 
way of distribution of Cu(i) 's is opposite for each unitary 
model. Thus the total Cum of valence bands is zero, but 
C sc = 2. In the case of A„ > X%, Cum's of each band 
vanishes, and C sc is also zero. For finite Ar, these two 
unitary models are hybridized and become a symplec- 
tic model. C sc of this hybridized system is quantized as 
C sc = 2 in the domain of the QSH state, while it vanishes 
in the domain of the usual SHI. With a finite magnetic 
field h x , similar hybridization occurs, but this breaks TR 
symmetry, and the model belongs to the topologically- 
trivial unitary class where Cum of each band is zero. In- 
terestingly, in the clean limit, C sc of this unitary model 
is still quantized, i.e, C sc = 2, while each pair of edge 
states opens a gap due to the hybridization by h x . (Note 
that C sc is well-defined even in TR breaking systems.) 

The localization length Xm{W,E) of a quasi one- 
dimensional tube of M-site circumference is calculated 
at energy E and disorder strength W by the transfer 
matrix method |3cj |. The M-dependence of the renor- 
malized localization length Ajw(W, E) = Xm(W, E}/M 
determines the localization/delocalization properties of 
the wavef unctions at energy E. 

Figures da-l)-(c-3) display A M (W,E) up to M = 24 
for several values of W as functions of E for disordered 
QSH system at h x = (symplectic ensemble). The case 
Xr = has been studied in the context of quantized 
anomalous Hall effect [2|| , and the two isolated extended 
states are identified by the M-independent Am(W, E). 
There are two energies at which extended states show 
up. As they merge together the extended states disap- 
pear. This is consistent with the scenario of levitation 
and pair annihilation of the two Cum 's of opposite signs. 
With finite A_r, the isolated extended states turn into fi- 
nite energy region of extended states, the width of which 
increases with Xr. Note that these two regions of ex- 
tended states approach as the disorder W increases. The 
gap in the density of states disappears already for W 
larger than ~ 3, but still the two regions are separated. 
When W is further increased, these two energy regions 
of extended states merge into one region, and eventually 
disappear. We have also checked that, in the ribbon ge- 
ometry, there appear extended gapless edge states even 
when there are no extended bulk states in the middle en- 
ergy region. Based on these results, we draw in Fig.^d) 
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FIG. 1: Renormalized localization length A M (W,E) = X M (W,E)/M for QSH with Ar = (a-l, 2,3), A R = 0.1 (b-1,2,3) and 
Ar = 0.2 (c-1, 2,3). The disorder strength W is increased from left to right as W = 5.0 (a,b,c-l), W = 7.0 (a,b,c-2) and W = 8.0 
(a,b,c-3). The other parameters are fixed as Aso = 0.3, A„ = 0.5 and h x — 0. (d) A localization/delocalization phase diagram 
obtained in the plane of energy E and disorder strength W . The red curve is the energy of the isolated extended states for 
Ar = 0, while the green and blue curves are the boundary of the energy region of the extended states for Ar = 0.1 and Ar = 0.2, 
respectively. 



a phase diagram depicting the location of extended states 
in the E-W plane. The red curve for \r = represents 
the trajectory of the isolated extended states in the uni- 
tary (QHE) case, while the phase boundaries between the 
localized and extended states are given by green curves 
for Ar = 0.1 and by blue ones for Ar — 0.2, respectively. 
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FIG. 2: Renormalized localization length Km(W,E) — 
Xm(W,E)/M for the unitary model with h x = 0.25 at (a) 
W = 3.0, (b) W = 4.0, (c) W = 5.0. The other parameters 
are Aso = 0.3, Ar = 0.1 and Ai, = 0.5. All the states are 
already localized since the U(l) Chern numbers for up and 
down spin bands cancel out. 

Let us confront some other cases with those discussed 
in Fig. 2] Figures GI a )-(c) display the curves Am(W,E) 
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FIG. 3: Renormalized localization length Am(W,E) = 
Xm(W, E)/M for SHI with A SO = 0.05, Ar, = 0.3 and A„ = 1.0 
at (a) W = 2.6, (b) W = 3.0 and (c) W = 3.4. 



for the unitary model with h x = 0.25 at (a) W = 3.0, 
(b) W = 4.0 and (c) W — 5.0. In this case, the extended 
states have already disappeared since Cum is zero for 
each of the split bands. Therefore the model is reduced 
to the trivial unitary class, where all states in two di- 
mensions are localized with any finite amount of disorder 
[25|. (It is also confirmed that the gapful edge states in 
the ribbon geometry are localized.) However, it should 
be noted that C sc for this system is quantized as C sc — 2 
in the clean limit. This means that, in TR breaking sys- 
tems, finite C sc does not protect extended states and that 
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protection would be closely related to the Kramers the- 
orem. One might think that the difference is due solely 
to symmetry, i.e., unitary v.s symplectic classes, and not 
to the topological property of the QSH state. It is there- 
fore important to compare with the simple symplectic 
model which belongs to the trivial Z2 classification. In 
our model, the SHI corresponds to this case, and Fig. |3| 
shows A M (W,E) for A so = 0.05, X R = 0.3, A„ = 1.0 at 
(a) W = 2.6, (b) W = 3.0 and (c) W = 3.4, respectively. 
It is evident here that the extended states disappear with 
much weaker disorder strength, and the two energy re- 
gions of extended states disappear without merging into 
a single one. 

The above two cases, i.e., Figs. [21 and strongly sug- 
gest that the localization behavior of QSH system in 
Fig. ^ is deeply influenced by the nontrivial topologi- 
cal aspect, and is distinguished from that of the usual 
symplectic class. In order to substantiate this expecta- 
tion, we have studied the critical property of the local- 
ization/delocalization transition of the QSH system. 

1 W = 5.0, E= 1.32, v= l.6l 
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FIG. 4: Scaling plot of the renormalized localization length 
Am(W,E) (see raw data inFig.HIb-1)) as a function of \E — 
E c \M 1/v for different energies E and M < 24 and for W = 5. 
The critical energy and exponent are estimated as E c — 1.32± 
0.005, v = 1.61 ± 0.10, respectively. 

Figure 01 summarizes the scaling analysis by displaying 
A M (W,E) = f(\E - E c \M l l u ) at W = 5.0 with v being 
the critical exponent for the divergence of the localization 
length. Data for various E and M (up to M = 24) are 
included and their collapse on a single curve indicates 
a reasonable one parameter scaling behavior, which si- 
multaneously determines v = 1.61 ± 0.10. We have also 
studied transition at higher disorder W = 7.0, and found 
v = 1.61 ± 0.10. This exponent should be confronted 
with that of the standard symplectic universality class 
v = 2.73 j2||, and that of the unitary model at strong 
magnetic field (with finite Cym) v = 2.33 [22]. The re- 
sult v = 1.6 obtained here is clearly distinct from both of 
these values, and suggests a new universality class for the 
symplectic ensemble with a non-trivial topological struc- 



ture. Intuitively, it is anticipated that the localization 
problem should be influenced by the finite spin Chern 
number. The construction of an effective field theory for 
this class is left for future investigations. 

In summary, we have studied the localiza- 
tion/delocalization problem of the QSH system (a 
symplectic ensemble with non-trivial topological prop- 
erty) by the transfer matrix method. The phase diagram 
shows levitation and pair annihilation of the two energy 
regions of extended states, analogous to that of the 
unitary model with finite Chern number (the integer 
quantum Hall effect). The critical exponent v for the 
localization/delocalization transition is estimated as 
v = 1.6, which is distinct from that of the standard 
symplectic universality class (v = 2.73) and that of the 
unitary class with non-zero Chern number [y = 2.33). 
This strongly suggests that the QSH system belongs to 
a new universality class characterized by a topological 
index such as the spin Chern number matrix or the Z2 
index. 
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